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The interplay between local, repulsive interactions and disorder acting only on one spin orienta¬ 
tion of lattice fermions (’’spin-dependent disorder”) is investigated. The nonmagnetic disorder vs. 
interaction phase diagram is computed using Dynamical Mean-Field Theory in combination with 
the geometric average over disorder. The latter determines the typical local density of states and 
is therefore sensitive to Anderson localization. The effect of spin-dependent disorder is found to be 
very different from that of conventional disorder. In particular, it destabilizes the metallic solution 
and leads to a novel spin-selective, localized phase at weak interactions and strong disorder. 
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I. INTRODUCTION 

Experiments with optically controlled cold atoms are 
able to turn new concepts and theoretical models into 
reality^ For example, the Bose-Hubbard Hamiltonian^ 
was realized and the superfluid-insulator transition was 
observed^ a fermionic Mott insulator was createdj^ An¬ 
derson localization was detected^ the spin Hall effect was 
measured)^ and exotic features such as fractional statis¬ 
tics were perceivedi^ In particular, such experiments 
make it possible to simulate models and discover phe¬ 
nomena which are absent in solid state physics. One such 
model, namely fermions on a lattice with spin-dependent 
hopping, was already realizedIn another experiment a 
spin-dependent periodic potential was formedi^ A possi¬ 
ble next step is to implement a lattice system with spin- 
dependent disorder by combining spin-dependent and lo¬ 
cal random potentials'^ 

Spin-dependent disorder may, in principle, be real¬ 
ized experimentally by focusing light beams with differ¬ 
ent polarization on an optical lattice, after having been 
scattered from a diffusive plate This leads to a spin- 
dependent, speckle-type effective random potential which 
is characterized by a short correlation length and pro¬ 
nounced statistical independence^^— and which acts dif¬ 
ferently on particles with different spin orientation in the 
ground state. Thereby spin-dependent speckle-type dis¬ 
order modifies the local one-particle energies, hopping 
amplitudes, and interparticle interaction potentials. 

In this paper we discuss the Anderson and Mott metal- 
insulator transitions (MIT) within the fermionic Hub¬ 
bard model in the presence of spin-dependent disorder. 
Surprisingly, up to now the influence of such a type of 
disorder on correlated lattice fermions has not yet been 
investigated in detail. Scalettar and collaboratorsii ana¬ 
lyzed the pairing instability in the attractive case, while 
the present authors studied the thermodynamical proper¬ 
ties in the repulsive case;^ By employing the dynamical 
mean-field theory (DMFT) with arithmetic average over 
the disorder it was showni^ that spin-dependent disor¬ 


der with a symmetric probability distribution function 
(PDF), Va{e) = T’ct(— e), yields a finite magnetization in 
the repulsive Hubbard systems on a bipartite lattice away 
from half-filling. At half-filling this novel type of disor¬ 
der breaks the particle-hole and spin symmetries, but no 
magnetization appears. The absence of magnetization is 
explained by the fact that the symmetry of the spin re¬ 
solved local density of states (LDOS) Pa{oj) = Pa{—oj) is 
preserved in the interacting system, so that the number 
of states below the Fermi energy w = 0 is unchanged. 

There remains the questions how spin-dependent dis¬ 
order will affect the Mott-Hubbard transition and Ander¬ 
son localization, respectively. In particular, in the case 
of non-interacting fermions, particles with different spin 
are independent and Anderson localization can therefore 
occur only in that spin subsystem on which the disorder 
acts. In the following we will shown that when the in¬ 
teraction is turned on, metallic and Anderson insulating 
states for different spin subsystems coexist in a large part 
of the phase diagram. 

For this purpose we employ the DMFT together with 
the geometric average over the disorder, since the lat¬ 
ter is sensitive to Anderson localization even within a 
one-particle description— Employing the geometric 
rather than the arithmetic average corresponds to deter¬ 
mining the typical LDOSi^ within the DMFT. This ap¬ 
proach was recently used to calculate the paramagnetioii 
and antiferromagnetioi^ phase diagrams of the disor¬ 
dered Hubbard model as well as the phase diagram of 
the disordered Falicov-Kimball model— Thereby it was 
possible to determine the MIT due to disorder (Ander¬ 
son localization) and interactions (Mott-Hubbard transi¬ 
tion), respectively, as well as the rich transition scenario 
caused by the simultaneous presence of interactions and 
disorder, within a unified framework. A recent exper¬ 
iment on correlated, disordered fermionic cold atoms^S 
confirmed that the Anderson localization line indeed in¬ 
creases linearly with the interaction strength as predicted 

in Ref. [ 13 - 

In the following we will refer to disorder acting equally 
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on both spin subsystems as ’’conventional disorder" to 
distinguish it from the spin-dependent disorder discussed 
here, where different random one-particle potentials act 
on the individual spin subsystems. In this study we 
do not consider long-range antiferromagnetic order since 
current experiments are still at temperatures above the 
Neel temperaturei^ 

II. MODEL AND METHOD 

The system under consideration is modeled by the 
Anderson-Hubbard Hamiltonian at half filling with spin- 
dependent diagonal disorder 

H — ^ ^ A ^ ^ T 

ijcr i(7 

U > (1) 

where Qia (al^) is the annihilation (creation) operator of 
a fermion at site i, a = ±1/2 ='(, ^ is the z component 
of the spin, riia = is the particle number operator, 

tij is the hopping amplitude between site i and j, and U 
is the on-site repulsion. 

The local spin-dependent potentials Ua represent un¬ 
correlated random variables drawn from a PDF Va{x). 
In this study the spin-dependent disorder is modeled by 
a rectangular (box) PDF given by 

V.(-) = (I - |x|) , (2) 

where 0(y) is the Heaviside step function and A is the 
strength of the disorder. It means that particles with spin 
a =d move on a lattice with a random potential whereas 
particle with spin a =), propagate on a uniform lattice. In 
the absence of disorder the system has SU(2) spin sym¬ 
metry and, at half-filling, is also particle-hole symmet¬ 
ric. The spin-dependent disorder breaks both symmetries 
since the on-site energy tia is equivalently represented by 
a random local chemical potential fii = — (ei-|-±ei 4,)/2 and 
a random local Zeeman magnetic field hi = (ei-|-— eij,)/2^i^ 

The Anderson-Hubbard Hamiltonian © is solved 
within the DMFT.— Here we map the Hamiltonian ([T]) 
onto an ensemble of single-impurity Anderson models 

^IsiAM = ^ Co-rio. ± U+ 

(7 

(3) 

k(T kcr 

with random spin-dependent atomic energies drawn 
from the same PDF Vai^a) as in ©. The hybridization 
matrix elements V^cr and the dispersion relation Aka- of 
the bath fermions Cko- define the hybridization function 
r]a{uj) = X^k |rk(Tp/(w — E]^a) and are determined self- 
consistently in the following way: For each we calcu¬ 
late the impurity Green function Go-( 02 , 60 .), correspond¬ 
ing to dSl), and the LDOS po(o 2 ,eo) = —ImGo(o;, eo)/ 7 r. 


From this we obtain the geometrically averaged LDOS 
Pa(uj) = exp[(lnpo(o2,eo))], where (Q) = f deP^(e)Q(e) 
denotes the arithmetic average of Q(e). The local Green 
function, averaged over the ensemble, is determined by 
the Hilbert transform Gcr(uj) = f duj'pa{uj')/{uj — ui'). 
The local self-energy So.(o;) is then obtained form the k- 
integrated Dyson equation So-( 02 ) = 02 —77o.(o2) —1/Go.(o2), 
where describes the mean coupling of the sin¬ 

gle lattice site to the rest of the system in the DMFT. 
The DMFT equations are closed by the Hilbert trans¬ 
form Go(a 2 ) = / fi^7Vo(0/[a2 - ^ - ^ 0 ( 02 )], where fVo(^) 
is the noninteracting density of states (DOS), which re¬ 
lates the local Green function for a given lattice to the 
self-energy from the impurity model. These equations 
must be solved iteratively to reach self-consistency. For 
comparison we also calculate the arithmetically averaged 
LDOS pf^^{uj) = {p^iuj,e,)). 

We choose a model DOS, Ao(^) = 2^/0^ — 
where W = 2D is the bandwidth, and IF = 1 sets 
the energy unit. For this DOS the local Green func¬ 
tion and the hybridization function are simply related 
by Ga(uj) = D^rycr(w)/4.— Here, the system is assumed 
to be paramagnetic down to the zero temperature. The 
phase diagram and dynamical quantities for the model 
o are computed by solving the DMFT equations at zero 
temperature with the numerical renormalization group 
method (NRG) ^^'^^ in each iteration step. For this pur¬ 
pose we implemented the open source NRG Ljubljana 
codej^ 


III. RESULTS 


The main result of our investigation is the ground 
state phase diagram of the Anderson-Hubbard model 
with spin-dependent disorder shown in Fig. [TJ It dif¬ 
fers significantly from the corresponding phase diagram 
for conventional disorder Three different phase transi¬ 
tions take place: (i) a Mott-Hubbard type MIT for weak 
disorder A, where the correlation gap at the Fermi level 
opens in the LDOS for spin up and down, (ii) an An¬ 
derson MIT for weak interaction U, where the LDOS for 
spin up vanishes, and (iii) a Falicov-Kimball type MIT, 
where the DOS for spin down acquires a correlation gap. 
The correlated, disordered metallic phase borders on the 
Mott insulator and the Anderson localized phase, where 
the spin-up particles are localized while those with spin 
down are still itinerant. In the following we will refer to 
the latter phenomenon as ’’spin-selective localization". 
The novel spin-selective localized phase originates from 
the spin dependence of the disorder. In the following we 
describe the properties of these three phases and charac¬ 
terize the transitions between them in detail. 
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FIG. 1. Nonmagnetic ground state phase diagram of the 
Anderson-Hubbard model at half filling with spin-dependent 
disorder determined by DMFT with the typical local density 
of states (LDOS). The inset shows the phase diagram obtained 
from the arithmetically averaged LDOS. 


A. Metallic phase 

The correlated metallic phase is defined by a non-zero 
value of the geometrically averaged LDOS at the Fermi 
level, Pcr(O) 0, for both spin channels. In the metallic 
phase without disorder (A = 0) Pcr(O) corresponds to the 
non-interacting DOS iVo(O) for both spins and arbitrary 
U. In this case the Luttinger theorem is obeyed and Lan¬ 
dau quasiparticles at the Fermi level are well defined.— 
In the presence of spin-dependent disorder, Pcr(O) is re¬ 
duced for particles in both spin channels, but differently 
depending in the spin direction. Only for U = 0 are 
the spin down particles not affected by the disorder at 
all. The reduction of Pct(O) due to an increase of A or 
[7 is shown in Fig. [2j In contrast to the case of conven¬ 
tional disorder— where a sufficiently strong interaction 
was found to protect the quasiparticles from decaying 
due to impurity scattering, Pcr(O) now always decreases 
with U for any finite A. Spin-dependent disorder vio¬ 
lates the Luttinger pinning condition at any U, and the 
local interaction cannot restore the Landau quasiparticle 
picture. As we will discuss later, this originates from the 
SU(2) spin symmetry breaking by this type of disorder. 


B. Mott-Hubbard MIT and the coexistence regime 

In the Mott insulating phase the LDOS vanishes at 
the Fermi level, Pcr(O) = 0. This is seen in the upper 
panel in Fig. [21 which leads to Hubbard subbands for 
stronger interactions U, cf. Fig. [S] The MIT at weak 
disorder, A < 0.8, is accompanied with hysteresis (upper 
panel in Fig. [2]), and the transition lines {U) and 
A^^{U) are tilted to the left (Fig. [T|). These transition 



FIG. 2. Local density of states at the Fermi level as a function 
of interaction strength U (upper panel) and disorder A (lower 
panel). 


lines terminate at a single critical point close to C/ ~ 0.6 
and A Rs 0.8. In contrast to conventional disorder^ 
where the critical point is located at A ss 1.8, the coex¬ 
istence regime is found to be significantly smaller, and 
a crossover regime does not occur at all. These differ¬ 
ences must be attributed to the reduced spin-symmetry 
in the present problem. Namely, the quasiparticle cen¬ 
tral peak between the Hubbard subbands, which origi¬ 
nates from spin-flip scattering and thereby leads to the 
Kondo resonance, is destroyed when a Zeeman magnetic 
field is applied— The induced finite magnetization re¬ 
duces the number of scattering states at the Fermi level 
in different spin channels. Although in the present prob¬ 
lem the magnetization is zero^ the density of states in 
the two spin bands differ, cf. Fig. [S] Altogether, spin- 
dependent disorder reduces the metallicity, forcing the 
interacting system into an insulating state, in contrast to 
the spin-independent case. 

The spin asymmetry of the bands leads to an un¬ 
equal population of particles with different spin in ex¬ 
tended states. This is seen in Figs. 01171 where we il¬ 
lustrate the dependence of the density n„ of fermions 
with spins a (upper insets) and of the polarization 
p = |(n-|- — n:i)/(n'[ + ni)\ (lower insets) on the disorder 
strength A for different interactions U. Without inter¬ 
action only the spin-up band is affected by the disorder. 
Due to the interaction between up and down particles 
the effect of the disorder is transmitted to the spin down 
band. The difference in the occupation of extended states 
is reduced by the interaction t/, which leads to a decrease 
of the polarization. 
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FIG. 3. Local density of states for different values of the 
interaction U at disorder strength A = 0.5. Upper panel: 
disordered metallic phase away from the transition line 
middle panel: disordered metallic phase close to the transi¬ 
tion line A^**; lower panel: disordered Mott insulator at the 
critical interaction strength U — 0.6. 




FIG. 4. Upper panel: Density of particles with spin up and 
spin down, respectively, in extended states at U — 0.25 as a 
function of disorder strength A. Lower panel: Polarization of 
extended states at U — 0.25 as a function of disorder. 


C. Spin-selective Anderson localized phase 

In the weakly interacting case, U < 0.55, a spin- 
selective Anderson localization is found for A > A^(t7), 
as shown in Fig. [T] Interestingly, we find that only par¬ 
ticles with spin up become localized when A > A^(t7), 
cf. Fig. [TJ whereas the particles with spin down are still 
itinerant. The LDOS at the Fermi level is zero in the 
former case and non-zero in the latter, as shown in the 
lower panel of Fig. [5] by the blue and green lines, respec¬ 
tively. On physical grounds this spin-selective localized 



FIG. 5. Upper panel: Density of particles with spin up and 
spin down, respectively, in extended states at U = 0.5 as a 
function of disorder strength A. Lower panel: Polarization of 
extended states at [/ = 0.5 as a function of disorder. 



FIG. 6. Upper panel: Density of particles with spin up and 
spin down, respectively, in extended states at U — 0.75 as a 
function of disorder strength A. Lower panel: Polarization of 
extended states at U = 0.75 as a function of disorder. 


phase may be effectively interpreted within a Falicov- 
Kimball model^ where spinless lattice fermions interact 
with immobile particles. In the present case, the par¬ 
ticles with spin up are localized due to the disorder A 
and form an immobile subsystem on which the parti¬ 
cles with spin down are scattered due to the interaction 
U. This spin-selective localization implies the absence 
of spin-up quasiparticle states around the Fermi level. 
However, Hubbard subbands are formed in the LDOS for 
both spins at higher energies, as is shown in the middle 
panel of Fig. |51 Again we note that although the disorder 
acts only on spin-up particles, the interaction transfers 
the effect of the disorder also to the spin-down particles. 
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FIG. 7. Upper panel: Density of particles with spin up and 
spin down, respectively, in extended states at U = 1.0 as a 
function of disorder strength A. Lower panel: Polarization of 
extended states at U — 1.0 as a function of disorder. 


This is evident form Figs. [5] and [51 where it is seen that 
some states with spin down become localized. Indeed, 
by setting the hybridization function for particles with 
spin up to zero, one can show that the DMFT action 
for the Hubbard model reduces to that for the Falicov- 
Kimball model^. This justifies our interpretation of the 
spin-selective Anderson localized phase in terms of the 
Falicov-Kimball model. 


D. Disordered Mott insulator 

For strong disorder, A > 0.8, we find a transition 
from the spin-selective Anderson localized phase to a 
disordered Mott insulator upon increasing U (’’Falicov- 
Kimball type MIT”). Namely, there is no hysteresis and 
the spin-down band is completely split off for Uc > 0.55. 
Above this value the LDOS at the Fermi level for spin 
down particles vanishes, cf. Fig. [5] The transition line 
(U) is seen to be vertical for disorder A > 0.8. 
This means that the localized spin-up fermions now play 
the role of the immobile particles in the Falicov-Kimball 
model. By further increasing U for a fixed A a Mott 
gap proportional to U develops, cf. the lower panel in 
Fig. [51 On the insulating side the integral over the LDOS 
changes with U although A is held constant. This is 
due to the Hubbard type coupling between spin up and 
down particles, leading to changes in both hybridization 
functions Therefore, we conclude that although 

the MIT above A w 0.8 is of the Falicov-Kimball type, 
the phase on the right hand side of the phase boundary 
in Fig. |T] is still a disordered Mott insulator. In other 
words, there is no phase boundary between weakly and 


strongly disordered Mott insulators. 



FIG. 8. LDOS at disorder strength A = 1 for different values 
of the interaction U. Upper panel: disordered metallic phase; 
middle panel: transition point from the disordered metallic 
to the spin-selective Anderson localized phase; lower panel: 
disordered Mott insnlator. 


IV. CONCLUSIONS AND OUTLOOK 

In this paper we discussed the properties of corre¬ 
lated lattice fermions in the presence of disorder acting 
only on spin-up particles. We computed the paramag¬ 
netic ground state phase diagram and one-particle quan- 
tites such as the spin resolved spectral functions and the 
particle densities per spin in the extended states. Fur¬ 
thermore, we identified three different phase transitions 
which separate four distinct phases: a correlated metal, 
a disordered Mott insulator, a spin-selective localized 
phase, and a coexistence regime. In the spin-selective lo¬ 
calized phase at strong disorder and weak interaction the 
spin-up particles are localized, whereas spin-down parti¬ 
cles are itinerant. 

The phase diagram and the properties of such a cor¬ 
related disordered system can be studied experimentally 
when a polarization dependent disorder is applied to an 
optical lattice filled with fermionic cold atoms. 

In the future it will be interesting to investigate also 
antiferromagnetic long range order, the effect of temper¬ 
ature, and the effect of a gradual turning-on of the spin 
asymmetry of the disorder. Work along these lines is in 
progress. 
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